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FOREWORD 

The  research  reported  in  this  report  was  performed  by  Willow  Run  Laboratories, 
a  unit  of  The  University  of  Michigan's  Institute  of  Science  and  Technology.  The  work 
was  supported  by  the  Advanced  Research  Projects  Agency  of  the  Department  of  De¬ 
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the  current  Project  Officer  is  Frank  Ruskey.  The  Willow  Run  Laboratories'  num¬ 
ber  for  this  report  is  10018-7-P. 
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period  include  Ben  Drake,  Philip  Jackson.  Ralph  Mitchel,  Roger  Turpening.  Robert 
Vincent,  and  Thomas  Wagner.  Roger  Turpening  is  Principal  Investigator,  while 
Philip  Jackson  Is  Principal  Scientist.  Section  3  of  this  report  was  specifically  con¬ 
tributed  by  Robert  Vincent.  Section  4  by  Ralph  Mitchel,  and  Section  5.2  by  Ben 
Drake.  The  remainder  of  the  report  was  written  by  Philip  Jackson. 

The  views  and  conclusions  contained  in  this  document  are  those  of  the  authors 
and  should  not  be  Interpreted  ns  necessarily  representing  the  official  policies, 
either  expressed  or  implied,  of  the  Advanced  Research  Projects  Agency  or  the  U.  S. 
Government. 
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ABSTRACT 

An  investigation  designed  to  utilize  simultaneously  aerial  photographs,  micro¬ 
wave  radar,  and  multispectral  scanners  for  geologic  studies  was  commenced.  In 
conjunction  with  ground-truth  investigations,  we  used  aerial  photographs  and  data 
acquired  from  both  passive  multispectral  scanning  and  multispectral,  multipolar, 
microwave  radar  systems  in  the  attempt  to  construct  geologic  and  lithological  maps 
applicable  to  the  selection  of  tunnel  sites.  Specific  areas  discussed  in  this  report 
include:  (1)  Vincent's  image- ratioing  technique,  which  provides  for  the  determina¬ 
tion  of  Si02  and  iron  oxide  content  in  specific  geologic  formations,  thus  allowing 
discrimination  between  different  types  of  surface  and  subsurface  features:  (2)  a 
discussion  of  different  types  of  coherent  optical  processing  techniques  available 
for  the  enhancement  of  geologic  structures  in  radar  imagery:  (3)  the  potential  for 
using  both  high- resolution  L-band  and  X-band  like-  and  cross-polarized  synthetic- 
aperture  radar  imagery  in  geologic  investigations;  and  (4)  a  discussion  of  the 
present  state  of  the  investigation,  with  recommendations  for  the  improvement  of 
present  techniques  and  the  implementation  of  additional  diagnostic  aids. 
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TUNNEL-SITE  SEI.ECTION  BY  REMOTE  SENSING  TECHNIQUES 
Semiannual  Report 
1  April  Through  30  September  1971 

1 

SUMMARY 

The  purpose  of  this  contract  is  to  evaluate  signal-enhancement,  data-reduction  techniques 
for  an  airborne  system  which  would  utilize  simultaneously  aerial  photographs,  microwave  radar, 
and  niultispectral  scanning  systems  (the  latter  operating  in  the  0.32-  to  13.5-spectral  region). 
The  data  acquired  from  such  a  system  are  to  be  utilized  in  the  construction  of  both  surface  and 
subsurface  geologlp  maps  which  could  be  used  to  predict  subsurface  conditions  in  possible  tun¬ 
neling  sites.  Ground-truth  investigations  are  to  be  conducted  to  verify  our  analyses  of  the  re¬ 
mote  sensing  data. 

For  these  evaluationa  we  are  using  both  presently  existing  and  newly  acquired  remote 
sensing  Imagery  collected  by  Willow  Run  Laboratories.  We  are  presently  improving  and  re¬ 
fining  our  techniques  for  collecting  and  analyzing  remote  sensing  data.  Recent  advancements 
in  the  area  Include:  (1)  the  development  of  Vincent's  image- ratloing  method,  which  utilizes 
multlspectral  Imagery  to  detect  differences  in  the  SiC^  and  iron  oxide  content  of  geologic  for¬ 
mations;  and  (2)  development  of  high- resolution  L-band  radar,  which  allows  us  to  acquire  llke- 
and  cross-polarized  Images  of  both  X-  and  L-band.  These  developments  in  remote  sensing 
techniques  appear  to  be  applicable  to  the  identification  of  geologic  features.  Such  new  techniques 
may,  in  fact,  be  required  before  the  use  of  remote  sensing  in  tunnel-site  selection  is  feasible. 

We  performed  a  literature  search  upon  the  application  of  airborne  radar  and  niultispectral 
scanning  systems  to  geologic  studies  in  order  to  avoid  duplication  of  effort,  to  build  upon  work 
performed  elsewhere,  and  to  acquaint  ourselves  generally  with  the  problems  involved  in  using 
remote  sensing  for  geologic  purposes.  We  have  now  acquainted  ourselves  with  much  of  the  sig¬ 
nificant  work  in  the  geologic  applications  of  remote  sensing  and  with  most  sources  of  current 
contributions. 

At  Willow  Run  Laboratories,  compositional  remote  sensing  (the  determination  of  differences 
in  lithology)  has  been  advanced  by  two  new  techniques  of  image  ratioing  performed  upon  multi- 
spectral  data.  These  techniques  allow  the  discrimination  between  silicate  and  nonsilicate  geo¬ 
logic  formations  and  the  determination  of  the  amount  of  ferric  oxide  in  given  geologic  features. 
Reflectance  and  emlttance  calibrations  are  now  being  made  on  rock  samples  to  refine  and  quan¬ 
tify  these  methods  further. 

Several  coherent  optical  proce*  dng  techniques  were  applied  to  available  radar  imagery  in 

r 

an  attempt  to  reduce  the  amount  oi  otal  information  and  enhance  those  features  applicable  to 
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the  identification  and  interpretation  of  geologic  features.  The  most  successful  technique  in¬ 
volves  an  image  enhancement  method  which  removes  the  low  spatial  frequency  portion  of  the 
image's  spectrum  bv  means  of  a  spatial  filter.  This  method  resulted  in  the  capability  to  delin¬ 
eate  more  easily  the  structural  edges  and  lines  separating  differently- reflecting  terrain  regions. 

During  this  period,  we  have  also  interpreted  high-resolution  L-band  radar  imagery  along 
with  high- resolution  X-band  imagery.  With  like-  and  cross-polarized  returns  in  each  band,  we 
have  the  capability  for  both  dual-spectral  and  dual-polarity  returns.  Interpretation  of  these 
radar  images  from  one  varied  geologic  region  revealed  ^itferences  in  the  Imagery  related  to 
the  texture  of  the  surface  and  to  the  substratum  under  a  thin  surface  layer.  iVj  differences 
shown  in  the  L-band  and  X-band  imagery  designate  those  regions  for  which  ground-truth  inves¬ 
tigations  should  be  initiated. 

With  the  use  of  compositional  remote  sensing  and  the  addition  of  multlspectral  and  multi¬ 
point-  radar,  we  presently  anticipate  that  remote  sensing  will  be  an  economically  useful  tech¬ 
nique  to  guide  field  investigations  of  potential  tunnel  sites.  We  now  need  to  quantity  more 
accurately  our  application  of  compositional  remote  sensing,  to  relate  multispectral  scanning, 
radar,  and  aerial  photography  more  completely  to  the  compositional  technique,  and  to  investigate 
computer  selection  and  color  coding  for  the  problem  of  tunnel-site  selection. 
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2 

INTRODUCTION 

Sections  3  through  5  of  this  report  discuss  specific  tasks  performed  during  this  reporting 
period.  Sections  3  and  4  present  our  progress  in  the  area  of  Imagery  Optimization  Analysis 
(Section  1.3  of  the  contract  work  statement),  while  Section  5  consists  of  a  Feasibility  Analysis 
(Section  1.4  of  the  work  statement).  Section  6  includes  an  evaluation  of  the  work  already  per¬ 
formed,  a  discussion  of  its  potential  for  geologic  interpretation,  and  a  summary  of  our  plans 
for  continued  investigations  into  the  application  of  remote  sensing  techniques  to  the  identification 
of  geologic  structures.  Based  upon  the  increase  in  knowledge  of  the  geologic  features  of  an  area 
acquired  through  remote  sensing  techniques,  we  hope  to  improve  the  present  methods  at  our 
disposal  in  order  to  aid  economically  in  the  determination  of  the  most  advantageous  sites  for 
tunnels. 
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3 

IMAGERY  OPTIMIZATION  ANALYSIS 

Gt'n logic  remote  sensing  has  at  its  disposal  two  new  techniques  of  image  ratioing  which 
are  applicable  (1)  to  the  discrimination  between  silicate  and  nonsilicate  geologic  formations 
and  (2)  to  the  determination  of  the  amount  of  ferric  oxide  in  given  geologic  features.  The  first 
technique  allows  the  discrimination  of  silicate  rocks  on  the  basis  of  a  ratio  of  the  radiances 
in  two  infrared  (1R)  channels.  Utilizing  25  rock  samples  selected  from  R.  J.  P.  Lyon  |1 1,  Fig. 

1  shows  the  ratio  between  channel  1  (8.2  to  10.9  pm)  and  channel  2  (9.4  to  12.1  pm)  and  the 
amount  of  Si02  in  percentage  of  weight,  with  the  ratio  symbolized  by  Rj  2 ■  The  crude  correla¬ 
tion  between  R(  2  -*d  the  percentage  of  Si02  is  evidenced  by  the  straight  curve  on  the  graph. 

The  wavelengths  of  the  silicon-oxygen  reststrahlen  bands  are  monitored  by  Rj  2.  The  more 
felsic  rocks  (those  with  greater  SiC>2  content),  generally  speaking,  have  reststrahlen  bands 
(relatively  large  departures  from  unit  emisstvity)  at  shorter  wavelengths  than  do  the  mafic 
rocks  (those  with  lesser  SiO,,  content),  and  therefore  Rj  2  increases  with  decreasing  Si02- 

The  second  image -ratioing  technique  involves  a  ratio  of  a  channel  in  the  visible  green  re¬ 
gion  of  the  spectrum  (channel  5,  0.50  to  0.52  pm)  and  one  in  the  reflective  1R  (channel  7,  0.74 
to  0.85  /  m).  This  ratio  can  show  the  rise  in  reflectance  of  ferric  oxides  from  the  green  to  the 

reflective  lit  wavelengths.  This  spectral  phenomenon  is  caused  by  electronic  transitions  of 
3  * 

the  Fe  ion  in  the  ferric  oxides  hematite  and  limonite. 

The  first  technique,  utilizing  the  ratio  of  two  1R  channels.  Is  demonstrated  in  Figs.  2-5, 
which  display  infrared  images  of  specially  processed  scanner  data  gathered  by  The  University 
of  Michigan  aircraft  at  an  altitude  of  3000  ft  over  a  sard  quarry  at  Mill  Creek.  Oklahoma,  and 
over  an  area  near  Pisgah  Crater  In  Southern  California  (also  at  3000-ft  altitude).  The  Mill 
Creek  data  were  gathered  on  June  25,  1970,  at  1000  hours  local  time,  while  the  data  from  Pisgah 
Crater  were  collected  at  0800  hours  local  time  on  October  30,  1970.  Each  of  the  figures  includes 
analog  infrared  images  of  two  single-wavelength  channels  of  data,  plus  an  image  of  the  ratio  of 
the  radiances  in  the  two  channels.  The  ratio  image  reflects  changes  in  the  chemical  composition 
of  the  rock  targets,  primarily  changes  in  their  Si02  content. 

Figure  2  shows  images  of  the  Mill  Creek  sand  quarry.  In  the  ratio  image  (Fig.  2c),  all 
features  that  appear  dark  (except  for  wavy  noise  lines)  are  either  exposed  quartz-sand  or  sand¬ 
stone.  The  outcroppings  around  the  lake  in  the  lower  part  of  the  images  are  sandstone.  Geolo¬ 
gists  from  both  the  United  States  Geological  Survey  and  The  University  of  Michigan  have  found 
almost  |ierfect  correlation  between  the  dark  regions  indicated  on  the  ratio  Image  and  quartz- 
sand  >r  sandstone  This  type  of  discrimination  should  prove  useful  in  explorations  for  construc¬ 
tion  materials,  such  as  sand  and  gravel. 
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FIGURE  1.  RATIO  BETWEEN  RADIANCES  IN  CHANNEL  1  AND  CHANNEL  2 
VERSUS  THE  PERCENTAGE  OF  WEIGHT  OF  Si02  IN  25  ROCK  SAMPLES. 

Samples  taken  from  Ref.  [  1| . 
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(a)  Channel  1:  8.2-10.9/im 


<b)  Channel  2:  9.4-12.1  /im 


tr)  Ratio  linage  of  Channel  1  over  2 

NOT  REPRODUCIBLE 

nu'  FIGURE  2.  DISCRIMINATION  OF  ACIDIC  SILICATES 

NEAR  MILL  CREEK.  OKLAHOMA  (SAND  QUARRY) 
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Figure  3  Is  '.ne  first  third  of  a  north-to-south  flight  over  Pisgah  Crater,  located  in  the 
right-hand  portion  of  eacli  of  the  three  images.  The  right  half  of  these  images  consists  primar¬ 
ily  of  basaltic  lava  flows  which  erupted  from  Pisgah  Crater  in  the  late  Pleistocene  or  Recent 
epochs.  These  relatively  young  lava  flows  (approximately  10,000  to  100,000  years  old)  erupted 
in  three  phases.  Although  similar  In  chemical  composition,  they  have  quite  different  surface 
textures,  varying  from  ropy  to  blocky.  The  older  langlomerate  in  the  left  half  of  the  images 
consists  primarily  of  rock  fragments  from  both  the  basaltic  lava  and  the  felsic  mountains  on 
both  sides  of  this  region.  North  of  the  highway,  parts  of  the  second  eruptive  phase  of  the  Pisgah 
lava  are  partially  covered  by  alluvium. 

With  the  exception  of  Pisgah  Crater,  topographically  the  highest  feature  in  the  images, 
the  single-channel  images  show  very  little  temperature  variation  across  the  scene.  The  ratio 
Image  (Fig.  3c),  however,  shows  considerable  variation  In  the  wavelengths  of  the  reststrahlen 
bands  over  the  same  region.  Lighter  and  darker  areas  on  the  ratio  image  Indicate  reststrahlen 
bands  at  longer  and  shorter  wavelengths,  respectively.  The  mafic  rocks,  which  contain  lower 
amounts  of  SlOg.  should  appear  bright  (Indicating  a  high  ratio),  and  the  felsic  rocks,  containing 
a  higher  SiOg  content,  should  appear  dark  (evidence  of  a  low  ratio).  In  the  ratio  image  (Fig.  3c), 
the  basaltic  lava  appears  bright  In  contrast  to  the  darker  tones  of  the  fanglomerate  and  gravel. 

The  basaltic  peninsula  near  the  center  of  Fig.  3c  can  easily  be  discriminated  from  the  surrounding 
fanglomerate  and  gravel,  whereas  the  single-channel  images  (Fig.  3a  and  3b)  show  very  low 
thermal  contrast  for  this  feature.  In  addition,  the  wind-blown  sand  extending  over  the  lava  just 
south  of  the  peninsula  shows  up  readily  In  the  ratio  image,  but  is  not  distinguishable  in  the  sin¬ 
gle-channel  Images.  In  the  ratio  Image,  the  three  bright  spots  on  Pisgah  Crater  are  the  only 
experimental  artifacts:  these  resulted  from  the  much  higher  temperature  on  the  eievated  sun¬ 
ward  slopes,  which  caused  clipping  of  the  signal  In  one  of  the  channels.  Patches  of  exposed 
basaltic  lava  north  of  the  highway  are  also  evident  in  Fig.  3c. 

Figure  4  shows  the  westernmost  half  of  a  west-to-east  flight  line  over  a  region  south-south¬ 
east  of  Pisgah  Crater.  From  left  to  right,  the  warm  region  (the  light  area  in  the  single-channel 
images)  in  the  extreme  western  part  is  a  dacitic  mountain:  the  adjacent  colder  region  is  fanglom¬ 
erate  and  gravel  (grading  eastwardly  to  smaller  fragments  and  sand):  and  the  broad,  warm  re¬ 
gion  is  the  basaltic  Sunshine  lava  flow.  Just  east  of  the  lava  flow  is  alluvium,  followed  by  playa 
deposits  (primarily  clay  and  carbonates)  and  the  southern  tip  of  the  basaltic  Pisgah  lava  flow. 

Once  again,  the  ratio  Image  (Fig.  4c)  shows  emissivity  variations  Indicative  of  rock  type.  The 
felsic  mountains  and  fanglomerate  appear  darker  than  the  mafic  lava,  and  the  playa  material 
is  contrasted  sharply  against  the  alluvium.  Because  the  mountains  are  warmer  than  other  parts 
of  the  scene,  computational  temperature  corrections  should  heighten  the  contrast  between  dacite 
and  basalt  by  lowering  the  ratio  illustrated  in  the  ratio  image.  In  other  words.  Fig.  4c  shows 
a  higher  ratio  for  this  section  of  the  image  than  would  be  expected  on  the  basis  of  the  wavelengths 
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(c)  Ratio  Image  of  Channel  1  over  2 


FIGURE  3  ANALOG  INFRARED  IMAGES  OF  FLIGHT  LINE  1.  SECTION  A  NEAR  PISGAH 

CRATER.  CALIFORNIA 
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♦ 

N 


(a)  Channel  1:  8.2-10.9 


(b)  Channel  2:  9.4-12.1  /im 

NOT  REPRODUCIBLE 


(c)  Ratio  Image  of  Channel  I  over  2 


FIGURE  4.  ANALOG  INFRARED  IMAGES  OF  FLIGHT  LINE  2.  SECTION  A.  NEAR  PISGAH 

CRATER.  CALIFORNIA 
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of  the  reststrahlen  bands  alone.  In  contrast,  the  patchy  appearance  of  the  playa  is  caused  solely 
by  emissivitv  variations,  since  the  temperature  is  uniform  across  the  floor  of  Lavic  Lake,  of 
which  this  section  is  a  part. 

Figure  5.  a  continuation  of  Fig.  4,  shows  the  western  half  of  the  flight  line  across  Lavic 
Lake.  From  west  to  east,  the  first  apparent  feature  is  the  Pisgah  basaltic  lava  with  playa 
material  below  it.  then  alluvium,  a  basaltic  mountain,  and  an  intermittent  stream  bed.  now  dry, 
which  has  eroded  down  through  the  surrounding  mountains  of  basalt.  Once  again,  on  the  ratio 
image  (Fig.  5c).  the  southern  end  of  the  Pisgah  lava  flow  appears  similar  to  the  playa  material, 
and  the  playa-ailuvium  boundary  is  sharply  demarcated.  On  the  northern  side  of  the  basaltic 
mountain,  the  small,  dark  regions  near  the  center  of  the  image  have  been  confirmed  by  recent 
field  trip  to  be  felsic  outcroppings.  The  two  most  prominent  dark  areas  in  the  stream  bed  are 
regions  of  feisic  volcanic  tuff,  which  contain  minor  amounts  of  malachite,  a  copper  carbonate. 
The  basaitic  mountain  in  the  lower  right  comer  of  Fig.  5c  appears  bright,  as  did  the  basaltic 
lava  flow  in  Fig.  4c.  As  in  the  previous  ratio  images,  in  Fig.  5c  the  felsic  rocks  and  rock 
fragments  comprising  the  alluvium  appear  to  have  lower  ratios  (and  therefore  appear  darker 
in  the  ratio  image)  than  do  the  mafic  rocks. 

These  images  show  that  by  ratioing  images  in  the  two  IR  channels  8.2  to  10.9  pm  and  9.4 
to  12.1  |im,  one  can  enhance  emissivity  variations  in  the  presence  of  temperature  gradients 
on  the  order  of  10°C  or  15°C  without  computational  temperature  corrections.  Even  in  shadowed 
regions,  compositional  information  is  still  present  in  the  ratio  images.  Although  temperature- 
corrected  ratio  maps  are  available  in  digital  form,  they  do  not  show  the  contrasts  in  geologic 
formations  as  dramatically  as  the  analog  infrared  ratio  images  do. 

Figures  6  and  7  illustrate  the  second  ratioing  method  —  i.e.,  the  ratioing  of  channels  in 
the  visible  green  and  the  reflective  IR  spectral  regions.  Figure  6  shows  the  same  area  as  that 
of  Fig.  4  (south-southeast  of  Pisgah  Crater),  except  that  the  data  have  been  recorded  in  channel 
5  (0.50  to  0.52  pm)  and  channel  7  (0.74  to  0.85  pm).  The  ratio  of  these  two  spectral  regions 
is  sensitive  to  the  presence  of  iron  oxides,  with  the  darker  areas  corresponding  to  greater 
contents  of  iron  oxide.  In  the  ratio  image  of  Fig.  6,  the  andesite  dikes  in  the  dacite  porphyry 
mountain  are  very  clear'v  outlined,  while  they  are  hardly  discernible  in  the  single-channel  maps 
of  the  area.  The  greater  amount  of  iron  oxide  in  the  basaltic  lava  (as  contrasted  to  that  in  the 
aiiuvium)  is  also  cleariy  illustrated.  Figure  7,  which  covers  the  same  area  as  that  of  Fig.  3, 
shows  that  the  lava  flows  around  Pisgah  Crater  vary  in  iron  oxide  content,  a  result  in  agree¬ 
ment  with  our  analysis  of  samples  recently  collected  in  the  field  from  each  of  the  three  erup¬ 
tive  phases. 

These  ratioing  techniques  have  just  recently  been  developed  and  are  presently  being  im¬ 
proved.  Under  the  present  contract,  new  laboratory  data  are  being  measured,  and  the  resuits 
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♦ 


(a)  Channel  1:  8.2-10.9  iim 


(b)  Channel  2:  9.4-12.1  fim 


NOT  REPRODUCIBLE 


(c)  Ratio  Image  of  Channel  1  over  2 


FIGURE  5.  ANALOG  INFRARED  IMAGES  OF  FLIGHT  LINE  2.  SECTION  B.  NEAR  P1SGAH 

CRATER.  CALIFORNIA 
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(a)  Channel  5:  0.50-0.52  11m 


(b)  Channel  7:  0.74-0.85 


(c)  Ratio  Image  of  Channel  5  over  7 


FIGURE  6.  ANALOG  VISIBLE  AND  NEAR-INFRARED  IMAGES  OF  FLIGHT  LINE  2.  SECTION 

A.  NEAR  PISGAH  CRATER.  CALIFORNIA 


NOT  REPRODUCIBLE 
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not  Reproducible 


(a)  Channel  5:  0.50-0.52  /im 


(b)  Channel  7:  0.74-0.85 /im 


(c)  Ratio  Image  of  Channel  5  over  7 


FIGURE  7.  ANALOG  VISIBLE  AND  NEAR-INFRARED  IMAGES  OF  FLIGHT  LINE  1.  SECTION 

A.  NEAR  PISGAH  CRATER.  CALIFORNIA 
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of  these  measurements  should  greatly  improve  our  ability  to  relate  the  measured  Rj  ^  *°  the 
chemical  parameters  of  geologic  fargets.  Either  the  relationship  between  SiOg  and  Rj  ^  t*1 
Fig.  1  will  be  improved  or  a  different  chemical  or  mineralogical  parameter  better  related  to 
Rj  2  will  be  found.  In  regard  to  the  R^  ^  ratio  technique,  several  samples  from  the  Pisgah 
Crater  test  site  are  presently  undergoing  laboratory  analyses  (both  spectral  and  chemical)  in 
an  effort  to  prove  quantitatively  the  relationship  between  spectral  reflectivity  features  and  the 
iron  oxide  content  of  exposed  rock  f  .faces.  It  seems  likely  that  the  new  developments  in 
image- ratioing  techniques  will  lead  to  improved  methods  for  rapid  geologic  mapping  and  litho¬ 
logic  differentiation. 
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4 

OPTICAL  PROCESSING  FOR  ENHANCEMENT  OF  GEOLOGIC  FEATURES 

Coherent  radar  imagery  has  shown  considerable  usefulness  for  all-weather  mapping  of 
terrain.  The  complex  scattering  mechanism  in  radar  imaging  is  strongly  affected  by  geomor- 
phic  factors  and  tends  to  enhance  subtle  topography  and  vegetation  [2| .  The  abiiity  of  a  radar 
image  to  reveai  smaii  tectonic  features  and  lineaments  has  been  recognized  by  others:  Rowan 
and  Cannon  [3|  have  utilized  radar  imagery  to  locate  both  previously  mapped  and  newiy  dis¬ 
covered  faults  in  the  Mill  Creek.  Oklahoma  area. 

Although  radar  imagery  effectively  displays  a  iarge  number  of  geomorphic  features,  it  stiii 
absorbs  a  coniderable  amount  of  a  geologist's  activity  to  locate  and  map  tectonic  features. 
Particularly  in  view  of  the  image  coverage  afforded  by  a  radar  sensor,  whereby  geoiogical  map¬ 
ping  caii  be  extended  over  much  of  the  earth's  surface,  the  present  capacity  to  coiiect  radar 
imagery  is  far  beyond  that  of  geologists  to  examine,  analyze,  and  interpret  It  adequately. 

This  problem  could  be  alleviated  by  processing  the  imagery  to  reduce  the  amount  to  be 
examined,  allowing  the  geologist  to  use  his  talents  for  Interpretation  of  the  structural  signa¬ 
tures.  Since  the  images  are  two-dimensional,  it  is  natural  to  consider  a  form  of  optical  pro¬ 
cessing  which  can  operate  upon  the  imagery  without  transformation. 

Coherent  optical  systems,  in  particular,  are  suitable  for  interpreting  radar  imagery. 

Along  with  the  capability  of  analyzing  each  point  of  the  image  simultaneously,  coherent  optical 
image -processing  systems  possess  a  physically  accessible  Fourier  transform  piane.  allowing 
the  spectrum  of  an  image  to  be  modified.  Thus,  linear  filtering  operations  can  be  performed 
upon  the  imagery  very  rapidly.  Examples  of  the  filtering  operations  commonly  performed  by 
a  coherent  optical -processing  system  are  directional  filtering,  matched  filtering,  and  image 
enhancement. 

A  series  of  experiments  was  performed  on  some  high-resolution,  synthetic -aperture  radar 
imagery  collected  by  The  University  of  Michigan  to  evaluate  the  usefulness  of  coherent  optical 
processing  as  an  aid  to  the  geologist.  The  area  mapped  was  near  Mill  Creek.  Oklahoma,  an 
area  underlain  by  gently  folded,  severely  faulted  pre -Pennsylvanian  rocks,  parts  of  which  are 
covered  by  Pennsylvanian  sandstone  and  conglomerates  of  the  lower  part  of  the  Pontotoc  Group. 
The  major  faults  and  folds  in  these  pre -Pennsylvanian  rocks  run  northwest -southeast.  In  many 
regions,  the  conglomerates  lack  the  extensive  structures  which  might  reveai  offsetting.  Rubble 
along  a  fault  trace  appears  similar  to  the  parent  rock,  making  the  recognition  of  a  fauit  diffi¬ 
cult,  even  with  meticulous  field  work. 
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Optical  directional  filtering  experiments  were  attempted  on  this  imagery  with  little  surress. 
Narrow  angular  wedges  placed  in  the  optical  system  to  modify  the  spectrum  of  the  radar  imag- 
erv  produced  p<x>r  discrimination.  Although  the  filtered  image  had  poor  contrast,  it  did  con¬ 
tain  many  streaks  identifying  lineaments.  However,  this  type  of  filtering  did  not  effectively 
reduce  he  information  to  be  interpreted.  Similar  results  were  obtained  when  a  two-dimension¬ 
al  Gaussian  slit  was  used  to  modify  the  spectrum  of  the  imagery. 

In  an  effort  to  obtain  better  discrimination,  optical  matched  filtering  was  attempted.  The 
mathematical  model  that  we  utilized  for  the  radar  energy  reflected  from  a  fault  trace  and  re 
corded  on  a  radar  image  was  as  follows: 


i(x.  y) 

where  i 
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the  radar  imagery  reflected  from  a  fault  trace  and  :•(  corded  on  a  radar  Image 
coordinate  across  width  of  fault  trace 
coordinate  along  length  of  fault  trace 
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width  of  the  fault  trace 
length  of  the  fault  trace 
path  taken  by  a  fault  trace 

center  of  fault  trace  as  a  function  of  the  location  along  Its  length 
center  of  fault  length 


2,\  and  x  are  written  as  functions  of  the  length  y.  since  the  trace  may  not  be  of  constant  width 
or  exactly  straight. 


We  defined  probability  measures  over  xq,  A,  and  a.  along  with  the  associated  probability 

density  functions  fx  (x)  and  f^(x)  over  x.  as  well  as  fa(y)  over  y.  XQ(y)  and  My)  are  random 

processes:  i.e..  for  each  y  they  are  random  variables.  We  assumed  that  fx  (x)  Is  symmetric 

around  some  value  x  (v)  M  (the  centroid  of  a  fault  trace  In  the  x-directlon);  that  f.  (x)  Is 
o  A 

symmetric  at  about  A(y)  A  (average  width  of  a  fault  trace);  and  that  f  (y)  is  symmetric  at 

a 

about  some  mean  value  a  a  (average  length  of  a  fault  trace). 

The  probability  i(x.  y)  IQ  is  the  probability  (P)  that  xq-A<x<xo  +  A  and  -a  -  y  ■:  a. 
which  ran  lie  computed  to  be: 

Pfi(x.  y  lo|  JffA(-x)  -  fA(x)|  .  Fx  (x)J[F_a(y)-  Fa(y)| 
jfF^(x)- FA(x)|  .  fx  <x)J[F_a(y)-  Fa(y)| 


where  F  the  probability  distribution  function. 
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A  particuiar  case  shouid.  however,  be  noted;  when  the  fault  trace  is  narrow,  having  a  con¬ 
stant  width  (2A  2f).  where  2f  is  on  the  order  of  the  radar  image's  resolution,  we  have: 

Pfifx.  y)  ij  fx  (x)(2c)[Fgy)  -  Fa(y)| 
o 

That  is,  the  probability  of  the  fault  trace's  existence  at  a  particuiar  point  is  defined  in  terms  of 
the  probability  density  function  of  the  random  variabie  defining  its  path. 

For  the  case  of  uniform  noise,  a  coherent  optical  correlator  using  a  matched  fiiter  in  the 
form  of  a  complex  spatial  fiiter  can  display  an  output  light  amplitude  distribution  proportional 
to  the  cross  correlation  between  a  reference  signal  (stored  on  the  fiiter)  and  the  input  image. 

The  results  of  the  experiment  utilizing  optical  matched  filtering  were  mixed.  Some  of  the  known 
faults  were  not  detected.  However,  with  the  fiiter  described  above,  we  were  ai)io  to  locate  the 
contact  separating  the  Deese  and  Dornick  Hiiis  formations. 

The  last  experiment  performed  was  image  enhancement,  designed  to  eliminate  from  the 
radar  image  that  information  which  is  unnecessary  to  geologic  interpretation.  The  iow  spatial 
frequency  portion  of  the  image's  spectrum  was  removed  by  means  of  a  spatiai  fiiter.  In  the 
processed  image,  structural  edges  and  iines  dividing  differently  reflecting  terrain  regions 
were  emphasized.  As  a  result,  this  experiment  allowed  geologic  interpretation  with  much  less 
time  and  effort  than  that  required  by  the  unprocessed  pictures. 
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5 

FEASIBILITY  ANALYSIS 

5.1.  DATA 

Multifrequency  ;ind  multipolar  radar  has  been  used  to  investigate  the  potential  for  identifica¬ 
tion  and  interpretation  of  geologic  features  in  terms  of  roughness,  general  texture,  composition, 
and  the  presence  of  features  on  the  order  of  a  wavelength  directly  below  the  surface.  The  re¬ 
sulting  imagery  is  now  being  used  as  an  aid  in  geoiogic  interpretation. 

Under  a  NASA  contract,  during  the  iast  week  in  August  The  University  of  Michigan  flew 
both  X-hand  and  L-band  radar  over  the  Pisgah  Crater  region  in  California.  Eight  signals, 
iike-  and  cross -polarization  of  both  initial  vertical  and  horizontal  orientation,  were  recorded 
for  each  frequency  on  this  flight.  Although  the  actual  resolution  of  the  eight  images  is  classi¬ 
fied.  the  images  of  the  Pisgah  Crater  did  reveal  differences  in  the  geologic  formations  of  the 
region.  A  dried  lake,  several  young  lava  flows  of  both  blocky  and  ropy  types,  alluvial  deposits, 
minor  exposures  of  intrusive  igneous  rocks,  and  some  apparent  drainage  features  were  de¬ 
picted  to  varying  extents  on  the  different  images.  Although  the  geologic  causes  of  some  of  the 
differences  between  different  images  can  be  inferred,  other  differences  cannot  be  explained 
from  examination  of  the  images  alone,  and  therefore  a  ground-truth  investigation  could  be  very 
helpful. 

Although  it  is  anticipated  that  the  images  can  be  presented  in  a  later  technical  report,  they 
cannot,  unfortunately,  be  reproduced  at  this  time.  However,  a  few  of  the  anomalies  in  the  de¬ 
gree  and  type  of  geologic  discrimination  in  different  images  can  be  pointed  out  as  examples. 

(1)  The  X-band  does  not  discriminate  between  lava  flows  and  a  thin  layer  of  iava  chips 
lying  on  silt,  while  the  L-band  does  discriminate  between  these  two  features. 

(2)  Differences  between  piaya  and  alluvium  show  up  more  in  X-band.  Possibly  this  means 
that  the  niiuvium  is  sluiiiowand  that  the  L-l)and  is  reflecting  from  a  slightly  buried 
feature  of  suitable  roughness. 

(3)  A  dendritic  area,  which  may  lie  a  drainage  feature,  shows  up  in  X-band.  but  is  greativ 
diminished  on  the  L-band  imagery.  This  pattern  does  not  show  up  in  muitispectral 
imagery,  but  the  reason  is  unknown. 

Perhaps  the  most  important  result  of  utilizing  radar  images  for  geologic  investigations  is 
that  they  have  provided  information  not  seen  previously  on  muitispectral  images  or  aerial 
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photographs,  although  ground-truth  information  should  bo  obtained  to  explain  specific  areas  lo¬ 
cated  on  the  radar  images.  In  the  comparison  of  the  multispectrai  and  radar  images  of  the 
Pisgah  Crater  region,  a  total  of  32  questions  were  raised  concerning  the  meaning  of  the  differ¬ 
ences  indiscrimination  of  various  geologic  features  on  specific  images. 

5.2.  EFFICIENT  UTILIZATION  OF  REMOTE  SENSING  IMAGERY 

Inspection  of  radar  imagery  of  various  polarizations  and  wavelengths  indicates  that  the 
maximum  amount  of  lithological,  geomorphological.  and  structural  information  is  obtained  when 
both  like-  and  cross -polarized  images  of  two  distinctly  different  wavelengths  are  available. 

!  ^ss  information  can  be  acquired  from  both  like-  and  cross -polarized  images  of  only  one  wave¬ 
length.  When  one  has  either  like-  or  cross -polarized  images  of  one  wavelength  (but  not  both), 
a  minimum  amount  of  information  is  obtainable,  although  even  in  this  case  significant  geologic 
data  can  be  derived. 

Very  little.  If  any,  lithological,  geomorphological.  or  structural  Information  can  be  ac¬ 
quired  from  aerial  photographs,  multispectrai  imagery,  or  radar  imagery  of  densely  built-up 
urban  areas,  unless  surface  geomorphological  patterns  (such  as  a  string  of  lakes  or  a  linear 
valley)  indicative  of  subsurface  lithology  can  be  seen.  However,  infrequently  limited  amounts 
of  such  information  can  be  obtained  from  aerial  photographs,  multispectrai  imagery,  or  radar 
imagery  of  cultivated  fields. 

Each  of  the  three  remote  sensors  (aerial  photography,  multispectrai  scanners,  and  radar) 
commonly  yields  lithological,  geomorphological.  and  structural  data  that  are  not  present  on  the 
Imagery  of  the  other  two  sensors.  If  at  all  possible,  aerial  photographs,  both  black-and-white 
and  color,  and  multispectrai  and  radar  imagery  of  the  same  area  should  be  obtained. 

5.3.  RADAR  RESOLUTION  FOR  DIFFERENT  WAVELENGTHS  AND  EVALUATION  OF  X-BAND 
IMAGERY  BEFORE  DEGRADING 

Regardless  of  the  resolution  of  the  Image,  the  amount  of  radar  reflection  is  dependent  upon 
the  relationship  between  the  wavelength  and  the  surface -roughness  scale  of  the  reflecting  sur¬ 
face.  For  example,  road  gravel  would  diffusely  reflect  1-cm  radar  waves  and  specularly  re¬ 
flect  30-cm  waves.  If  Images  were  degraded  to  30-  by  30-ft  resolution,  a  large  area  of  g*  avel 
would  appear  light  (relatively  high  reflection)  at  1  cm  and  dark  (relatively  little  reflection)  at 
30  cm.  Both  experimental  and  theoretical  analyses  show  that  the  different  wavelengths  are 
diagnostically  useful  even  when  degraded. 

Nondegraded  X-  and  L-band  imagery  has  been  interpreted  for  geologic  purposes,  and 
tentatively,  it  appears  that  more  pertinent  geologic  information  can  be  obtained  from  those 
images  with  higher  resolution. 
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:>.4.  LITERATURE  SEARCH 

The  literature  search  was  performed  to  study  the  application  of  airborne  radar  and  mult i - 
spectral  scanning  systems  to  geologic  studies.  We  employed  sources  from  the  Department  of 
Defense  and  the  National  Aeronautics  and  Space  Administration;  bibliographies  were  purchased 
or  otherwise  acquired,  pertinent  current  journals  and  government  reports  scanned,  and  basic 
texts  studied.  The  pur|K>scs  of  the  search  were  to  avoid  duplication  of  effort,  to  build  upon  and 
include  in  our  evaluation  work  which  had  been  performed  elsewhere,  and  to  acquaint  ourselves 
with  the  problems  and  (Kissible  approaches  to  their  solutions. 

Although  we  have  acquired  a  general  knowledge  of  the  state  of  the  art  and  of  the  sources  of 
current  contributions,  we  still  have  much  to  learn.  For  example,  photogeology  has  been  de¬ 
veloped  to  a  highly  sophisticated  discipline  requiring  much  experience.  Since  our  project  is  con¬ 
cerned  with  developing  and  evaluating  an  extension  of  photogeology,  we  must  be  certain  to  use 
properly  the  techniques  and  potential  of  photogeology  for  extended  applications  with  multispec- 
tral  and  radar  data.  Photogeology  is  only  one  of  the  areas  in  which  continuous  literature  study 
is  required. 

A  secondary  objective  of  the  literature  search  has  been  to  study  the  geologic  problems  of 
tunneling.  Thus  far.  we  have  studied  the  general  problems  that  tunneling  presents  in  order  to 
tie  aware  Inith  of  the  geologic  features  to  lie  sought  and  those  to  be  avoided  when  tunneling  sites 
are  selected. 

In  subsequent  reporting  periods,  our  emphasis  will  be  on  current  literature  dealing  with  the 
general  area  of  remote  sensing  as  applied  to  geologic  research. 

a. 5.  T11K.  FKASiniLlTY  OF  EXTENDING  25-db  TONE  AND  TEXTURE  OF  RADAR  RECORD¬ 
INGS  TO  60-70  db 

Essentially,  the  signal  film  recorded  from  synthetic-aperture  radar  is  in  a  frequency  with 
relatively  wide  range  (  60  db)  (Drown  and  Porcello.  1969  [4 1;  Cutrona.  et  al.,  ’.966  [5 1).  When 
the  output  of  the  radar  optical  processor  is  viewed  directly,  the  wide  dynamic  range  is  retained: 

,i  very  bright  spot  in  the  output  will  stun  the  eye  just  as  shining  a  flashlight  diiectly  Into  it  will. 
When  an  output  film  is  made  of  the  radar  scene,  the  bright  s  pots  either  bloom,  or  the  exposure 
and  development  of  the  film  is  so  controlled  that  the  brightness  is  reduced  to  about  20  db. 

A  new  method  of  making  a  hologram  of  the  radar  output  has  recently  been  developed  at 
W  i  1  low  Run  Laboratories,  so  that  the  hologram  may  be  viewed  through  rather  simple  optics 
which  retain  the  wide  dynamic  range.  Although  this  new  type  of  hologram  provides  an  extremely 
striking  effect,  it  is  not  evident  at  this  time  that  it  will  aid  geologic  interpretation. 

Under  another  contract,  the  Radar  and  Optics  Division  has  designed  and  is  installing  a 
digitizing  system  which  uses  an  image  dissector.  This  system  will  digitize  images  in  the  output 
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plane  of  the  radar  optical  processor,  so  that  the  wide  dynamic  range  will  he  available  for 
digitizing.  Time  of  completion  is  presently  estimated  to  he  from  two  to  three  months.  We 
anticipate  digitizing  radar  images  of  regions  being  studied  and  using  digital  computation  to 
perform  image  subtraction  or  ratiolng. 
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6 

PRESENT  STATE  OF  THE  INVESTIGATION.  POTENTIAL 
APPLICATIONS.  AND  PLANS 

In  most  of  the  specific  t;isks  delineated  in  the  contract,  significant  advances  have  been 
accomplished  at  the  Willow  Hun  Laboratories  during  this  reporting  period.  One  of  the  most 
important  advances  is  compositional  remote  sensing  made  possible  by  Vincent's  ratioing 
method,  which  was  onlv  hypothetical  at  the  time  the  work  for  this  contract  was  proposed.  This 
ratioing  technique  has  produced  a  new  dimension  in  geologic  remote  sensing. 

Only  certain  types  of  surface  lithologies  can  now  be  identified,  but  the  ratioing  method  is 
potentially  able  to  differentiate  between  ail  rock  types  in  which  the  reflectance  or  emittnnce 
spectra  differ.  An  example  of  the  potential  of  this  method  is  the  discovery  of  andesite  dikes  by 
compos  it  ionallv  sensing  their  ferric  oxide  content  in  a  NASA  test  region  which  previously  had 
been  thoroughly  investigated  both  by  remote  sensing  and  ground  mapping. 

The  extension  of  our  high- resolution  radar  to  L-band,  with  the  resulting  capability  of  pro¬ 
ducing  multispectrai  and  multipolar  images,  has  enabled  us  to  locate  geologic  features  not  pre¬ 
viously  identified.  We  plan  to  utili/.e  a  ratioing  technique  upon  these  radar  images  similar  to 
that  used  by  Vincent  upon  multispectrai  imagery. 

Rowan,  et  ai.  |6|  of  the  U.  S.  Geological  Survey  have  delineated  fault  zones  in  predawn 
thermal-lR  imagery,  an  accomplishment  which  may  have  significant  applications  to  the  selec¬ 
tion  of  tunneling  sites.  The  identification  of  fault  zones  and  use  of  the  thermal  inertia  properties 
of  different  rock  types  should  extend  the  capabilities  of  geologic  remote  sensing. 

In  addition,  digital  and  analog  computations,  such  as  ratioing.  color  coding,  and  possibly 
optical  processing,  should  aid  the  geologist  in  interpreting  remote  sensing  data. 

Our  task  now  is  not  only  to  improve  these  individual  approaches,  but  to  coordinate  them  to 
learn  both  their  limitations  and  their  potentials  in  selecting  accurately  and  economically  the 
best  tunnel  sites.  To  make  this  evaluation,  we  are  concentrating  on  geologic  areas  for  which 
we  already  have  much  imagery,  one  of  which  is  the  Pisgah  Crater  region  in  California.  Al¬ 
though  it  is  not  a  region  through  which  one  would  expect  to  build  tunnels,  it  does  contain  varied 
geologic  features  which  are  challenging  to  identify  and  analyze  and  which  might  be  encountered 
in  some  potential  tunnel  sites. 

Our  plans  arc  to  continue  investigating  and,  if  possible,  improving  specific  remote  sensing 
techniques.  We  intend  to  continue  the  detailed  geologic  interpretation  of  remote  sensing  data 
and  to  evaluate  our  interpretations  through  ground-truth  investigations  of  those  areas  of  signifi- 
i  nit  geologic  interest  for  which  we  already  have  extensive  data.  At  present,  these  areas  include 
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the  Pisgah  Crater  region,  the  Niagara  escarpment  region  in  southern  Ontario,  and  the  Arbuckle 
Mountains  near  Mill  Creek,  Oklahoma.  Possibly  after  conference  with  the  Project  Officer, 
NASA  flights  can  be  requested  for  the  purposes  of  this  project. 
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